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We present a new determination of the distance (and age) of the Galactic globular 
cluster 47 Tucanae (NGC 104) based on the fit of its white dwarf (WD) cooling sequence 
with the empirical fiducial sequence of local WD with known trigonometric parallax, 
following the method described in Renzini et al. (1996). Both the cluster and the 
local WDs were imaged with HST+WFPC2 using the same instrument setup. We 
obtained an apparent distance modulus of (m — M)y = 13.27 ± 0.14 consistent with 
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previous ground-based determinations and shorter than that found using HIPPARCOS 
subdwarfs. Coupling our distance determination with a new measure of the apparent 
magnitude of the main sequence turnoff, based on our HST data, we derive an age of 
13 ± 2.5 Gyr. 
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1. Introduction 

Galactic globular clusters (GCs) are the oldest 
systems for which ages can be accurately mea- 
sured using the stellar evolutionary clock. As 
such, they constrain the age of the Galaxy, hence 
of the Universe, and by extension the choice of 
cosmological parameters. The only observable 
quantity that can be used to measure the ab- 
solute age of a GC is the absolute magnitude of 
the main sequence turnoff (My). Other observ- 
ables are also sensitive to age, e.g. the horizontal 
branch (HB) morphology or the turnoff (TO) col- 
ors, but their use is dangerous since in addition to 
age they are sensitive to other parameters such as 
mass loss and convective efficiency, respectively. 

The TO luminosity of a GC is derived from 
the observed magnitude of TO stars, correcting 
for interstellar extinction, and applying the dis- 
tance modulus. Simple calculations show that 
the age (t) error budget is indeed dominated by 
the uncertainty in GC distances, with the simple 
rule St/t = <miod (Renzini 1991), with e.g., an 
uncertainty of 0.25 mag in distance modulus im- 
plying a ~ 25% uncertainty in age. Correspond- 
ingly, efforts at improving the accuracy of GC 
age determinations have focused on the distance 
issue. 

Although the direct measurement of the par- 
allax of several (perhaps all) GCs is on the hori- 
zon (SIM, GAIA), for the time being we are still 
bound to use indirect methods, based on stan- 
dard candles. In most recent years local subd- 
warfs have been widely used as standard candles, 
especially after trigonometric parallax data from 
the HIPPARCOS mission have become publicly 
available. 

Indeed, HIPPARCOS has produced a large 
sample of local subdwarfs with accurate paral- 
laxes, and the general perception has been that a 
sizable reduction of GC ages using the subdwarf 
method (Reid 1997, 1998; Gratton et al, 1997) 
was due to the new set of subdwarf parallaxes. 
Actually, the situation is somewhat more intri- 



cate. For subdwarfs in common with previous, 
ground-based parallaxes the new HIPPARCOS 
measurements give distance moduli that on av- 
erage are just 0.1 mag larger (see Fig. 2 in Reid 

1997) . This corresponds to a reduction in age 
of only ~ 10%, while claimed increases in GC 
distance moduli were ~ 0.3 mag, and the corre- 
sponding reductions over previous age determi- 
nations were up to ~ 30% (Reid 1997; Gratton 
et al. 1997). Therefore, HIPPARCOS parallaxes 
alone account for only ~ 1/3 of the claimed re- 
ductions, much of the effect coming instead from 
a variety of cumulative effects, from including 
subdwarfs and evolved stars with relatively un- 
certain parallax measurements (Reid 1997), to 
adopting a different metallicity scale (Gratton et 
al. 1997). Systematic differences of 1-2 Gyr are 
also introduced by the use of different theoret- 
ical My° —age relations. Moreover, contrasting 
results for M 92 (Pont et al. 1998, Reid 1997 
and Carretta et al. 2000) show how much care 
is required to properly deal with the many ob- 
servational biases and reddening and metallicity 
scales. In fact, the fortuitous agreement of the 
distance moduli derived by Pont et al. and by 
Carretta et al. is a spurious result, given by adop- 
tion of different reddening values and metallici- 
ties for the field subdwarfs used in the fitting. In 
summary, HIPPARCOS GC distances and ages 
still remain matter of debate. 

Every method of distance determination presents 
advantages and disadvantages. For example, GC 
distances derived from the subdwarf method re- 
quire the metallicity of the cluster and of each 
subdwarf to be specified, which makes the method 
prone to possible systematic biases of the dwarfs 
vs. giants metallicity scales (see e.g., King et al. 

1998) . Thus, a mismatch of ~ 0.25 dex in the 
subdwarfs vs. GC metallicity scales results in a 
~ 0.25 mag error in distance and in a 25% error 
in age (Renzini 1991). For these reasons, we be- 
lieve that a variety of methods should be used to 
derive GC distances, as discrepancies among the 
results should help to identify possible systematic 
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biases. 

As an alternative to the subdwarf method, 
Renzini et al. (1996, hereafter Paper I) used 
the DA white dwarf cooling curve as a fidu- 
cial to measure the distance to one of the near- 
est GCs, NGC 6752, obtaining a true modu- 
lus of 13.05 mag and an age of 14.5 Gyr (in- 
cluding He diffusion) and 15.5 Gyr (not includ- 
ing He diffusion). In this paper, we use the 
white dwarf (WD) method to derive the distance 
modulus to 47 Tuc, a globular cluster that is 
roughly a factor often more metal rich than NGC 
6752 ([Fe/H]=-0.70±0.03; Carretta & Gratton 
1997). In extending our approach to this higher 
metallicity, we will be able to determine whether 
there is an age-metallicity relationship among 
the Galactic GCs. The similarity of the Color- 
Magnitude diagram (CMD) locus of 47 Tuc to 
that of bulge globulars NGC 6528 and NGC 6553 
(Ortolani et al. 1995) will constrain the ages of 
bulge clusters of near solar metallicity, as well as 
of the field population of the Galactic bulge. In 
principle, the large difference in metallicity be- 
tween NGC 6752 and 47 Tuc should be useful in 
giving an independent constraint on the relation- 
ship between the V magnitude of the horizontal 
branch and [Fe/H]. Unfortunately, the horizontal 
branch of NGC 6752 is nearly vertical, so its V 
magnitude is not so well defined. On the other 
hand, 47 Tuc has only one peculiar RR-Lyrae 
star, but does have a well defined red HB at 
V = 14.10 (Kaluzny et al. 1997). These two 
clusters were selected for the application of the 
WD method for being the least reddened ones 
among the nearest GCs. 

2. Observations and Data Analysis 

A field located 6'5 West of the center of 47 Tuc 
was observed during Cycle-5 with the HST-WFPC2 
camera through the filters F336W, F555W, and 
F814W. Short exposures with the same instru- 
ment setup + long F439W band frames were 
taken during Cycle-7. Finally, some extra F814W 
frames of the same field, acquired for astrometric 



purposes, where kindly provided by I. King. All 
in all, the total exposure time was 26300s, 13800, 
15880 and 21949 seconds in F336W, F439W, 
F555W and F814W, respectively. Most of the 
observations were cosmic-ray split. For the lo- 
cal calibrators, we re-measured the magnitudes 
of the four DA WDs used in Paper I. In order 
to achieve a better definition of the local WD 
sequence, four new local WDs (2 DA + 2 DB) 
were observed in the same filters. One more star 
(WD1647+591) was observed but not used be- 
cause its location in the CMD is too red for its 
mass (0.69 M ). We conclude that either it is 
a binary, or that there must be an error in the 
mass determination and/or in the parallax. 

All the images were de-biased, dark corrected 
and flatfielded through the standard HST pipeline. 
Following Silbermann et al. (1996), we have 
masked out the vignetted region, the saturated 
and bad pixels and columns using a vignetting 
mask created by P.B. Stetson, together with the 
appropriate data quality file for each frame. 

2.1. Cluster data 

The photometric reduction of the frames in 
the field of 47 Tuc was carried out using the 
DAOPHOT II/ALLFRAME package (Stetson 1987, 
1994). Preliminary photometry was performed 
on each single frame in order to obtain an approx- 
imate list of stars for each of them. The coordi- 
nates of the common stars were used for an ac- 
curate spatial match among the different frames. 
With the correct coordinate tranformations, we 
then obtained a single image, which was the me- 
dian of all the frames, regardless of the filter. 
In this way we removed all the cosmic rays and 
obtain the highest signal-to-noise image for star 
identification. We ran the DAOPHOT /FIND 
routine on the median image and performed PSF 
fitting photometry on it, in order to obtain the 
deepest list of stellar objects, free from spurious 
detections. Finally, this list was given as input 
to ALLFRAME, for the simultaneous profile fit- 
ting photometry of all the individual frames. The 
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PSF we used were the WFPC2 reference PSFs 
extracted by P.B. Stetson (1995) from a large set 
of uncrowded and unsaturated images. 

For a consistent comparison between the 47 Tuc 
and the local WD stars, all the measured in- 
strumental magnitudes were scaled to a Is ex- 
posure. Aperture corrections for the cluster data 
were measured from a few isolated stars in each 
chip and filter. They were applied to the cluster 
stars in order to convert the PSF magnitudes to 
the same zero point as the local WDs, measured 
with an aperture of 0''5 radius. The stars ob- 
served with each of the WF chips were then put 
together, after removing the small chip-to-chip 
offset (Dolphin 2000). 

It is now well known that the WFPC2 detector 
is affected by a charge transfer efficiency (CTE) 
problem: a loss of charge during the CCD trans- 
fer phase. The amount of the loss depends on the 
position within the chip, the background counts, 
and has worsened with time since the installa- 
tion of the WFPC2 camera. In the last few years 
various estimates of the size of this effect have 
been made (e.g., Whitmore & Heyer 1997; Stet- 
son 1998; Whitmore, Heyer & Casertano 1999) 
that was sometimes mis-interpreted as a "long 
vs short anomaly" (Casertano & Mutchler 1998). 
However only very recently Dolphin (2000) pro- 
vided a formula that takes into account all the de- 
pendencies of the charge loss, including its varia- 
tion with the frame background (and therefore 
with the exposure time) and its increase with 
the epoch of the observation. The CTE problem 
is critical for our method, since it relies mainly 
on the photometry of faint stars in low back- 
ground (the background being especially low in 
the F336W frames). For this reason, we have 
paid special attention in applying the CTE cor- 
rection (Dolphin 2000) and in veryfing that it 
gives sound results. The size of the correction 
for a cluster WD in F336W ranges from a few 
hundredth up to 0.3 magnitudes, depending on 
its magnitude and position in the chip. However, 
the brightest cluster WDs have been corrected 



by less than 0.05 mag, and therefore if a residual 
CTE problem were present, a large difference be- 
tween the slope of the cluster and the local WD 
sequence would be present. The fact that no such 
difference is apparent (Fig. |6|) makes us confident 
about the adopted correction. Finally, a further 
correction, of the order of a few hundredth of 
a magnitude, was applied to the F336W magni- 
tudes to account for the time dependence of the 
UV throughput (Baggett & Gonzaga 1998). 

The instrumental CMDs containing the stars 
from the three WF, in the four color planes, are 
shown in Fig. [l]. The photometry from the PC 
chip was not used because no WDs were de- 
tected in this field. In what follows we will use 
the notation 777336, 771439, 777555, w-814 to indicate 
the Is instrumental magnitudes in the F336W, 
F439W, F555W and F814 filters, respectively. 
Only the stars detected in > 50% of the frames 
in each filter were accepted, and a further se- 
lection based on the magnitude errors and the 
PSF fitting parameters was made. Fig. [I] shows 
only the stars (~ 2300) in common to the four 
diagrams, for this reason the main sequence in 
the (771555, 777,555 — 777814) plane is truncated at 
a magnitude brighter than the 777555 limit mag- 
nitude. A WD sequence is clearly visible in all 
four diagrams, although it is better defined in 

the (777336, 777 336 ~ 777555 ) and (777 3 36 , 777 33 6 - 777 814 ) 

planes. The main sequence and subgiant branch 
of the Small Magellanic Cloud old population are 
also clearly visible at an intermediate location 
between the 47 Tuc main sequence and the WD 
sequence. A total of 21 cluster WDs were iden- 
tified. The WD sequence appears very narrow 
especially in the two CMDs using the 777336 mag- 
nitudes, which, due to the larger color baseline, 
are also the most useful one to distinguish be- 
tween DA and DB dwarfs ( cf. Fig. |2|) . Therefore, 
we can conclude that only DA WDs have been 
sampled in this 47 Tuc field, and therefore only 
the local WDs of the DA type will be used as 
calibrators for the distance determination. 
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2.2. Local WD data 

The local WD calibrators were selected for 
having a trigonometric parallax with accuracy 
better than 10% and a spectroscopically deter- 
mined mass as close as possible to 0.53 M & , 
the predicted mass for the globular cluster WDs 
(Renzini & Fusi Pecci 1988; see also Fig. 1 in 
Greggio and Renzini 1990). 

Every local WD was imaged by each of the 
3 WF chips. Aperture photometry was per- 
formed within a radius of 0.5 arcsec on each 
image, the chip-to-chip offset, as determined 
by Dolphin (2000), was removed, and then the 
three magnitudes of each star were averaged. 
The instrumental magnitudes were scaled to Is 
exposure and corrected for the decreasing UV 
throughput (m336 only). The absolute instru- 
mental magnitudes were then obtained by means 
of the trigonometric parallaxes reported in Ta- 
ble 1. Three of our local calibrating WDs were 
observed by HIPPARCOS, and their parallaxes, 
although slightly different from the ground-based 
ones, do not show any systematic zero point shift. 
These stars are WD0644+375, WD1327+276, 
and WD1917-077, and the percent differences 
{^hip — t^va) between the parallaxes in the HIP- 
PARCOS and the Van Altena et al. (1991) cat- 
alog are: +4%, —9%, —12%, respectively. Note 
that the latter is a DB WD, so it was not used in 
the fit. In addition, Vauclair et al. (1997) showed 
that the comparison between the parallaxes the 
16 DA WDs in common between the two catalogs 
reveals no systematic difference. 

The Lutz-Kelker correction, calculated accord- 
ing to the formula given in Carretta et al. (2000), 
amounts to a few hundredth of a magnitudes for 
all the stars except one, the DB WD WD0002+729, 
for which it is ~ 0.3 magnitudes. However, this 
star, being of DB type, was not used for the dis- 
tance determination. Table 1 also shows that the 
masses M\vd of a few of the selected local WDs 
differ appreciably from the ideal value, 0.53 M Q . 
In order to obtain a constant-mass sequence ap- 



propriate for the match with the cluster WDs, 
the correction 5(mag) = 2.4(M\vd — 0.53) as 
applied to the magnitudes of all the local WDs 
(Wood 1995; Bergeron, Wesemael & Beauchamp 
1995). Moreover, the WD masses listed in Ta- 
ble 1 were all (except WD2326+049) originally 
derived adopting a WD mass-gravity relation ap- 
propriate for WDs without an hydrogen enve- 
lope. However, as extensively discussed in Bra- 
gaglia, Renzini &; Bergeron (1995), there is now 
general consensus that "evolutionary" values for 
the mass of the hydrogen envelope should be 
used, i.e., ~ 10 _4 Mq for the WD masses of in- 
terest here (Napiwotzki, Green & Saffer, 1999). 
This implies that WD masses in Table 1 were 
underestimated by ~ 0.04M Q (Bragaglia et al. 
1995), and therefore such masses have been in- 
creased by this amount in producing the final 
fiducial WD sequence used in the distance de- 
termination. The magnitude correction for each 
calibrating WD (except WD2326+049 that was 
obtained with the correct assumption on the en- 
velope thickness) was then taken as 5(mag) = 
2.4(M WD + 0.04 - 0.53), where M WD are the 
masses in Table 1. 

The resulting fiducial WD sequences are shown 
in Fig. |^, including the best fit straight line in 
each individual CMD. This plot also shows the 
size of the Lutz-Kelker (solid vertical lines) and 
of the mass corrections (dotted lines) applied 
to each star. Note that the dispersion is very 
small and that the one reliable DB dwarf is lo- 
cated well outside of the DA sequence only in the 
(^336; m 336 — ^555) plane. As mentioned above, 
no star in the 47 Tuc CMDs matches the location 
of the DB sequence. Figure || shows the same 
CMDs but only for the WDs used for the dis- 
tance determination, with the estimated errors 
being also indicated. 

In summary, the WD fiducial cooling sequence 
relies on the spectroscopically determined gravi- 
ties, a theoretical mass-gravity relation corrected 
for the finite mass of the hydrogen envelope, 
and on the adopted mass of the cluster WDs, 
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0.53 ± 0.02Mq. This latter assumption is exten- 
sively justified in Paper I, using a series of con- 
vergent stellar evolution arguments that we be- 
lieve are quite solid. Of course, the direct spec- 
troscopic measure of the mass of cluster WDs 
would be quite reassuring. VLT spectroscopic 
observations of a sample of four HST selected 
WDs in the globular cluster NGC 6397 have re- 
cently confirmed their nature of WDs of the DA 
variety (Moehler et al. 2000), although the low 
S/N of the relatively short exposure spectra of 
these very faint stars did not allow an accurate 
determination of the mass. 

2.3. Photometric zero point check 

Globular cluster distance determinations have 
to rely on standard candles, at least until the di- 
rect measure of trigonometric parallaxes will be 
feasible with SIM/GAIA. The results of the cur- 
rent methods of distance determinations largely 
rely on the consistency between the observations 
of the standard candles and their cluster counter- 
parts. Therefore, ideally, the two samples should 
be measured precisely in the same conditions. 
However, in most cases this approach is not vi- 
able. In the present work, although we observed 
both the local and the cluster WDs with the same 
instrument setup, this fact alone does not guar- 
antee the necessary homogeneity, given that we 
are comparing the photometry from very deep 
exposures of a sample of extremely faint stars 
in a crowded field (i.e. object for which PSF 
fitting photometry is needed) with the aperture 
photometry of a few very bright stars in short- 
exposure frames of almost empty fields. There- 
fore, in order to verify that no residual photo- 
metric zero point difference can affect our con- 
clusions, we report here the results of some tests 
we carried out to check the consistency of our 
measurements, i.e., the linearity of WFPC2 over 
a range of ~ 10 magnitudes. 

The first, most obvious test of the instrumen- 
tal magnitudes is their calibration to the Lan- 
dolt system (following Holtzmann et al. 1995 and 



Dolphin 2000), and the comparison of the mea- 
sured calibrated magnitudes with other photom- 
etry available in the literature. We emphasize 
again that we will use only instrumental magni- 
tudes to estimate the distance of 47 Tuc. This 
test on the calibrated CMD has the sole pur- 
pose of checking for the presence of some sys- 
tematic errors in the instrumental magnitudes, 
that of course would show up in the calibrated 
ones. Since no other photometry is available for 
our field of 47 Tuc, it is possible only to compare 
our CMD with the fiducial lines of the previous 
photometry. Fig. ||a shows the comparison be- 
tween our calibrated (V, V — I) diagram and the 
fiducial lines of the photometries by Kaluzny et 
al. (1997; solid) and Ortolani (private communi- 
cation; dashed). There appears to be excellent 
agreement between Ortolani's photometry and 
this work, while Kaluzny's data seem to be sys- 
tematically bluer by ~ 0.07 magnitudes. Unfor- 
tunately we have no means of establishing which 
one of the two sets of ground-based photometry 
is the most correct. Fig. ^|b shows the compari- 
son among our (V, B — V) CMD and the fiducial 
sequences by Kaluzny et al. (1997; solid) and 
Hesser et al. (1987, dotted). In this case there 
is excellent agreement. Figure || makes us confi- 
dent about the zero point of our 771,439, "1555 and 
171814 magnitudes. In particular, the above test 
is also a check of systematic errors in the aper- 
ture corrections, that would affect our distance 
determination. 

We did not calibrate the 171,336 magnitudes 
because the F336W filter bandpass differs sig- 
nificantly from the Johnson U ; and, therefore, 
the calibration is unreliable, especially for the 
hottest /coldest stars. However, in order to check 
the zero point of the 771336 instrumental mag- 
nitudes, we performed an artificial star test on 
all the F336W frames. We did this experi- 
ments on the F336W frames only because our 
distance mainly relies on the 711336 magnitudes 
(due to the fact that the WD sequence in this 
band is less steep, better defined and well sep- 
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arated from the SMC main sequence) although 
the F336W S/N for the faintest WDs was not 
excellent. 100 artificial stars were added, in the 
same position, on each single frame, with mag- 
nitudes in the range of the cluster WDs. Com- 
plete DAOPHOTII / ALLFRAME reduction was 
then carried out on these frames with the same 
algorithm used for the original images, but man- 
ually adding these artificial stars to the star list. 
The output magnitudes of the artificial stars 
were then compared with the input ones. Fig. [| 
shows this comparison: the systematic error in 
the output magnitudes is consistent with zero. 
For this reason we can exclude the possibility of 
a displacement of the cluster WD sequence in 
the (m 3 3 6 ,m 3 3 6 - mgu) or (m 3 3 6 ,m 3 3 6 - 777555) 
diagrams, due to the migration of the magni- 
tudes towards brighter values, as a consequence 
of crowding (Stetson & Harris, 1988). Of course, 
as all the artificial star tests, this uses model 
PSFs that could be different from the real ones. 
Still, the fact that the mean difference between 
input and output magnitude is consistent with 
zero gives a good indication that no significant 
migration is present in the 771336 photometry of 
the faint stars. Note that the number of stars in 
Fig. H is not related to the completeness of our 
photometry because the artificial stars were not 
identified by the same finding algorithm used for 
the original ones. 

3. The distance to 47 Tuc 

Once the cooling sequences for the cluster and 
the local WDs are defined, the distance modu- 
lus is the vertical shift that makes the cluster 
sequence overlap the local one. We adopt a red- 
dening E(B — V) = 0.055, which is the mean 
value of those reported by Zinn (1980), Reed et 
al. (1988) and the Stromgren catalogue (Hauck 
&; Mermilliod 1990). From this value the corre- 
spondent interstellar absorption A\ was derived 
for the three bands, according to Table 12 in 
Holtzmann et al. (1995), and subtracted them 
from our magnitudes. Then we shifted the clus- 



ter WDs to match the local ones. In doing so, 
we adopted for the WD sequence in each clus- 
ter CMD the same slope determined for the local 
calibrators, leaving only the zero-point (i.e., the 
cluster distance modulus) as the adjustable pa- 
rameter. 

Fig. H shows the match between the two WD 
sequences in the three CMDs (msu, 777336— msu), 
(777336, "7336 - rn 555 ) and (777555,777555 - tt7 8 i 4 ). 
The distance modulus required to match the se- 
quences in each panel is shown by the label. The 
quoted errors are the internal ones, including the 
magnitude uncertainties both in the local and 
cluster WDs, and the errors in the adopted par- 
allaxes. The first two panels allow a more accu- 
rate distance determination, either because the 
WD sequence is less steep (left panel), or because 
of the smaller dispersion around the fiducial se- 
quence (central panel). Note that the scale of 
both axes has been kept the same in the three 
panels in order to show the difference in the se- 
quence slopes. The arrow in the lower left corner 
of each panel shows the direction of the redden- 
ing vector, whose size has been increased by a 
factor of three in order to make it more visible. 
Depending on the difference in slope between 
the reddening vector and the WD sequence, a 
fixed uncertainty in the cluster reddening corre- 
sponds to a different error in the distance mod- 
ulus for each panel. In fact, with an uncertainty 
in the cluster reddening AE(B — V) = ±0.02, 
the corresponding A(t77 — M)o are =f0.03, =f0.01 
and =p0.07 in the left, central and right panels 
of Fig. [| respectively. The very small varia- 
tion of the distance modulus derived from the 

(777814, 777 336 ~ ^8u) and the (t77 336 , 777 336 - 777555) 

CMDs is due to the fact that in these planes the 
WD sequence happens to have a slope very sim- 
ilar to that of the reddening vector (see arrows 
in Fig. ||). Therefore a change in the adopted 
reddening shifts the sequence almost parallel to 
itself. Adding the above quoted systematic un- 
certainties to the formal errors of the fit, and 
averaging the three results weighted by the in- 
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verse of their errors, we obtain a mean distance 
modulus of (m - M) = 13.09 ± 0.13. 

A further source of systematic error in this 
value comes from the adopted mean mass of the 
cluster WDs: M WD = 0.53 ± 0.02. According to 
the relation quoted above 5(mag) = 2.4(M\vd — 
0.53) an uncertainty of 0.02 in mass implies a 
0.05 systematic uncertainty in the distance mod- 
ulus. Our best estimate for the distance modulus 
of 47 Tuc is therefore (m - M) = 13.09 ± 0.14. 
Combining this value with the adopted redden- 
ing E(B - V) = 0.055 ± 0.02 and the extinction 
law Ay = 3.2E(B — V) the apparent distance 
modulus is (m - M)y = 13.27 ± 0.14. 

In summary, a formula that explicitly shows 
the dependence of the apparent distance modulus 
upon the two main assumptions, namely cluster 
WD masses and reddening, can be expressed as 
follows: 

(m-M)v = !3.09+3.2xE(B-V)-2x[E(B-V) -0.055] 
-2.4 x (Mwd - 0.53) + 2.4 x (< Mwd > -0.594) 

where < Mwd > is the true mean mass of the 
sample of local WDs, and M WD the true mass 
of the cluster WDs, and 0.594 is the mean mass 
of the sampled local DA WDs given in Table 1, 
including the correction for the thick hydrogen 
envelope. 

For consistency with Paper I, in our analy- 
sis we adopted the trigonometric parallaxes from 
the Van Altena et al. (1991) catalog. A new 
version of the same catalog, including some more 
stars and a new weighting criterion for the var- 
ious measurements has been released more re- 
cently (Van Altena et al. 1995). Due to the 
new weighting system, the parallaxes of two of 
the local DA WDs, namely WD0644+729 and 
WD2126+734, increased by 6% and 9%, corre- 
sponding to a decrease of 0.12 and 0.18 magni- 
tudes in their distance modulus, respectively. As 
a consequence, had we adopted such parallaxes, 
the distance modulus for 47 Tuc would have de- 
creased by ~ 0.1 magnitude. 



Table 2 shows a summary of the previous de- 
terminations of the distance of 47 Tuc, using a 
variety of methods, together with the adopted 
reddening. For a direct comparison among the 
different distances in Table 2 we reported all to 
the same reddening E(B — V) = 0.055, and the 
result is plotted in Figure 0, with the correspond- 
ing errorbars. Our distance determination lies at 
the lower end of the distribution. It is consis- 
tent within the errors with most previous deter- 
minations, but it is significantly shorter than the 
values reported by Reid (1998) and Carretta et 
al. (2000) using the subdwarf method based on 
HIPPARCOS parallaxes. The cause of this dis- 
crepancy remains to be understood. 

However, were the correct modulus as high as 
13.69, i.e., ~ 0.42 magnitudes larger than the 
value indicated by the WD method, then the 
inplied mass of the cluster WDs would be as 
low as 0.53 - 0.4/2.4 = O.36M , much less than 
indicated by stellar evolution theory. Alterna- 
tively, the mass of the local WDs (via Balmer 
line fitting) should have been underestimated by 
~ 0.17Mq, on average, which also looks quite 
implausible. 

Still assuming that the WD method is to 
blame, the discrepancy could arise from a mis- 
match between the cluster and the local WD 
magnitude scales. As shown by Fig. |6| the dis- 
tance moduli obtained from the three CMDs dif- 
fer appreciably one from another. This can only 
arise from a systematic error in the photometry 
of at least one of the three bands. Indeed, an er- 
ror of about Am555 = +0.05 would be sufficient 
to bring the three distance moduli into coinci- 
dence, at the value of (m — M)o = 13.10 found 
from the (obviously unvaried) (msi4, mzz&— m%i±) 
CMD. Similarly, an error of about Amsu = 
—0.08 would make the three distances coincide 
with the value of (m — M)q = 12.95 found from 
the middle panel of Fig. |(| Due to the still 
poorly known CTE correction, the presence of 
such systematic errors for the very faint cluster 
WDs cannot be ruled out. However, this would 
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not help reducing the discrepancy with the sub- 
dwarf distances given by Reid (1998) and Car- 
retta et al. (2000). Instead, a systematic er- 
ror of +0.2 in the 777,336 magnitudes only would 
bring the distance moduli derived from the two 
leftmost panels in Fig. || into coincidence with 
the value of (m — M)q = 13.34 found from 
the (771555,771555 — 777814) CMD. This would also 
solve the discrepancy with the Carretta et al. 
(2000) result. However, it would be more likely 
that there is a systematic error of the order of 
0.05 — 0.08 in the V and I instrumental magni- 
tudes, where the WDs are relatively fainter and 
the field is more crowded, rather than an error 
as large as 0.2 mag in the U band. 

Finally, it is worth noting that were a residual 
systematic error the same in all the magnitudes, 
then it would have no effect on the colors, and 
therefore it would just shift the derived distance 
moduli by the same amount (equal to the error 
itself) in all the CMDs. 

4. The age of 47 Tuc 

The present determination of the distance 
modulus, combined with the knowledge of the 
apparent magnitude of the TO and the global 
metallicity, allows us to derive the age of the clus- 
ter, via the turnoff luminosity My°. Figure || 
shows the cluster TO region in the CMD con- 
structed using three different color baselines^. 
The horizontal line shows our best estimate of the 
TO magnitude V TO = 17.65, while the two dot- 
ted lines define its error range estimate: ±0.10 
magnitude. Coupled with our distance deter- 
mination, this gives a TO absolute magnitude 
M$° = 17.65 - 13.27 = 4.38 ± 0.17. 

There is little scatter among recent determi- 
nation of the metallicity of 47 Tuc, and we adopt 

11 Note that, as mentioned above in the text, the calibrated 
U magnitude is not reliable at extreme colors because of 
the large difference between the F336W and the Johnson 
U bandpass. We illustrate it here because the relatively 
narrow width of the main sequence in this plane more 
clearly indicates the V magnitude of the turnoff point. 



[Fe/H]=-0.70 (Carretta & Gratton 1997). It is 
now generally believed that for metal poor halo 
stars and globular clusters the a-elements en- 
hancement requires values of [a/Fe] between 0.4 
and 0.6. However, following Gratton, Quarta & 
Ortolani (1986) we adopt [a/Fejw 0.3 for the 
abundance ratios in 47 Tuc. 

Figure |9| shows the My° vs age relation, for 
the adopted metallicity [M/H]=-0.5, from the 
most recent models available in the literature, 
and the values we derive are listed in Table 3. 
The first column in Table 3 gives the reference of 
the adopted model, the second and third give its 
metallicity and helium abundance, the fourth in- 
dicates whether the helium and heavy element 
diffusion was included in the models, and the 
next one gives the age corresponding to 47 Tuc. 
The average age obtained from the models that 
include atomic diffusion is 12.9 Gyr, while the 
others give 13.5. As expected (e.g., Castellani 
et al. 1997) the two values differ by about 0.6 
Gyr. Note however that the 1 Gyr discrepancy 
between the age from Cassisi et al. (1999) and 
Straniero et al. (1997) is mainly due to the use 
of a different equation of state in the two sets of 
models (OPAL the former, and a modified ver- 
sion of the Straniero 1988 the latter; see Cassisi 
et al. 1999, their Table 1). The age error cor- 
responding to the uncertainty in My is ~ 2.2 
Gyr. 

The last column of Table 3 shows the ages 
derived for NGC 6752 assuming = 17 A - 

13.27. Here V TO = 17 A (Penny & Dickens 1986) 
is the same value adopted in Paper I, while the 
distance modulus {m — M)y = 13.27 comes from 
the value found in Paper I ((771 — M)y = 13.17) 
corrected for the fact that in Paper I a thin hy- 
drogen envelope was assumed for the disk WDs, 
and therefore a correction of 2.4 x 0.04 is needed 
(see Sec. 2.2). As Table 1 clearly shows, the 
two clusters turn out to be coeval within the er- 
rors. A forthcoming paper will be devoted to a 
more direct comparison between the 47 Tuc and 
NGC 6752 turnoff locations, matching the cluster 
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WD sequences . 

5. Conclusions 

The distance of the Galactic globular clus- 
ter 47 Tuc as been measured by comparing its 
WD cooling sequence to the fiducial sequence 
obtained from a sample of local WDs of known 
trigonometric parallax. We derive an apparent 
distance modulus of (m — M)y = 13.27 ± 0.14, 
or (m — M)o = 13.09 ± 0.14, implying a turnoff 
absolute magnitude of = 17.65 - 13.27 = 

4.38±0.17 and therefore an age of ~ 13±2.2 Gyr, 
if models including atomic diffusion are used. 

We demonstrate again the feasibility of the 
method using deep optical WFPC2 imaging, that 
permits very accurate measurement of magni- 
tudes even for stars as faint as the GC WDs. The 
low dispersion in the cluster cooling sequence 
together with the fact that the method is not 
affected by uncertainties in the metallicity, nor 
in the absolute calibration of the photometry 
(instrumental magnitude system is consistently 
used) permit us to measure the distance moduli 
of the closest clusters with an error of 0.14 mag, 
which includes systematic uncertainties. 

The derived distance modulus is ~ 0.42 mag 
shorter than an extreme value obtained with the 
subdwarf fitting method. Such "long" distance 
for 47 Tuc would imply an implausibly low mass 
for the cluster WDs (~ 0.36M Q ), or a very large 
systematic error in the masses of the local WDs 
as obtained by Balmer-line fitting. 

Nevertheless, an important requirement for 
the application of this method is the linearity 
of the detector at the two extremes of the dy- 
namical range. HST allows one to obtain very 
small statistical photometric errors down to the 
faintest stars, but problems like the CTE loss in 
the WFPC2 may affect the accuracy of the mea- 
sures. Therefore, the use of the appropriate cor- 
rection for this effect can be crucial for the result. 
We applied a state of the art determination of the 
correction, but identify in the CTE problem the 



main possible source of residual systematic bias. 

However, to reconcile the "long" distance to 
47 Tuc, the mismatch between the cluster and 
the local WD sequences should be as high as ~ 
0.36 mag in U-I, which again we consider highly 
implausible. 
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Table 1 

THE LOCAL CALIBRATING WHITE DWARFS 



WD 


7T ± Alt 


Ref. 


M ± AM 


Ref. 


WD type 


WD0002+729 


0.0291 ± 0.0047 


1 


0.600 ± 0.030 


2 


DB 


WD0644+375 


0.0626 ± 0.0018 


1 


0.655 ±0.020 


3 


DA 


WD1327-083 


0.0611 ±0.0028 


1 


0.502 ± 0.017 


4 


DA 


WD1917-077 


0.1010 ±0.0026 


1 


0.550 ±0.050 


5 


DB 


WD1935+327 


0.0561 ± 0.0029 


1 


0.512 ±0.013 


6 


DA 


WD2126+734 


0.0433 ± 0.0035 


1 


0.513 ±0.012 


6 


DA 


WD2326+049 


0.0725 ± 0.0048 


1 


0.690 ± 0.025 


7 


DA 


WD2341+322 


0.0559 ± 0.0017 


1 


0.494 ± 0.021 


6 


DA 



References: (1) Van Altena et al. (1991); (2) Beauchamp (1995); (3) 
Bergeron, Saffer & Liebert (1992); (4) Bragaglia et al. (1995); (5) Oswalt et 
al. (1991); (6) Bragaglia & Bergeron (2000); (7) Bergeron et al. (1995) 



Table 2 

47 TUC DISTANCE DETERMINATIONS 



Method (m-M)v E(B-V) Reference 



HB fitting 


13.40 ±0.20 


0.040 


Hesser et al. (1987) 


Baade-Wesselink 


13.36 ±0.17 


0.040 


Storm et al. (1994) 


HB fitting 


13.45 ±0.07 


0.040 


Kaluzny et al. (1998) 


HB fitting 


13.50 ±0.05 


0.050 


Salaris & Weiss (1998) 


HB fitting 


13.43 ±0.15 


0.040 


Ferraro et al. (1999) 


RGB tip 


13.42 ±0.20 


0.040 


Ferraro et al. (2000) 


MS fitting 


13.69 ±0.15 


0.040 


Reid (1998) 


MS fitting 


13.55 ±0.09 


0.055 


Carretta et al (2000) 


WDs 


13.27 ±0.14 


0.055 


This work 
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Table 3 

THE AGE OF 47 TUC USING DIFFERENT MODELS 



Model [M/H] Y diffusion age (Gyr) NGC 6752 age 



Straniero et al. (1997) -0.60 0.230 yes 13.5 13.0 

Cassisi et al. (1999) -0.52 0.230 yes 12.4 11.7 

VandenBerg et al. (2000) -0.62 0.241 no 13.7 13.3 

Girardi et al. (2000) -0.50 0.245 no 13.3 13.5 

Salasnich et al. (2000) -0.50 0.245 yes 12.8 
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Fig. 1. — The instrumental CMD of the observed 47 Tuc field. Three main branches are clearly distin- 
guishable: the cluster MS on the right, the SMC MS and SGB in the middle, and the cluster white dwarf 
sequence on the left. 
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Fig. 2. — Instrumental CMDs of the local WDs used as calibrators. DA WDs are plotted as filled squares, 
DBs are the asterisks. A least square fit to the DA WDs is shown as a solid line. The vertical dotted 
lines show the size of the magnitude correction applied to the WDs whose masses differ from 0.53, while 
the solid vertical line is the size of the Lutz-Kelker correction. 
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Fig. 3. — Same as Fig. |2|, but only for the DA WDs and with errorbars. The latter include the internal 
photometric errors, as well as the uncertainties in the parallax and mass determinations listed in Table 1. 
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Fig. 4. — Comparison between our calibrated data and the fiducial lines of previously published pho- 
tometries. Left panel: V, V-I CMD compared with the fiducial lines from Kaluzny et al. (1997; solid) 
and Ortolani 2000 (dashed). Right panel: B, B-V CMD compared with Kaluzny et al. (1997; solid), and 
Hesser et al. (1987; dotted). 
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Fig. 5. — Artificial star test in the F336W band. The difference between the input and output magnitudes 
are plotted as a function of the input magnitude. 
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Fig. 6. — Match between the 21 WDs identified in 47 Tuc (open circles) and the local DA WDs (filled 
squares). The distance moduli required to match the two sequences in the three planes are shown in the 
labels, together with their formal error. The arrow on the lower left corner shows the direction of the 
reddening vector, whose size has been increased by a factor of three in order to make it more visible. The 
dotted line in the central and right panel shows the displacement of the cluster WD sequence if the mean 
distance modulus were adopted. 
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Fig. 7. — Determinations of the distance of 47 Tuc by different authors, in the same order as Table 2, 
from left to right. Different symbols refer to different methods: HB fitting (filled circles), the unique 
RR-Lyrae variable (open circle), RGB tip (filled square), subdwarf fitting (open squares) and this work 
(filled triangle). 
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Fig. 8. — The turnoff region in the calibrated CMD. 
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Fig. 9.— M$° vs age relation derived from different theoretical models. The horizontal dotted line shows 
the estimated location of My° for 47 Tuc, and the vertical ones the corresponding ages, according to the 
three models. 
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